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ABSTRACT 

We analyze the conditions for detection of CO(l-O) emission in the Large Magellanic Cloud (LMC), using 
the recently completed second NANTEN CO survey. In particular, we investigate correlations between CO 
integrated intensity and H I integrated intensity, peak brightness temperature, and line width at a resolution 
of 2.'6 (^40 pc). We find that significant H I column density (exceeding ~10^' cm"^) and peak brightness 
temperature (exceeding ~20 K) are necessary but not sufficient conditions for CO detection, with many regions 
of strong H I emission not associated with molecular clouds. The large scatter in CO intensities for a given 
H I intensity persists even when averaging on scales of >200 pc, indicating that the scatter is not solely due 
to local conversion of H I into H2 near GMCs. We focus on two possibilities to account for this scatter: either 
there exist spatial variations in the /(CO) to A^(H2) conversion factor, or a significant fraction of the atomic gas 
is not involved in molecular cloud formation. A weak tendency for CO emission to be suppressed for large H I 
linewidths supports the second hypothesis, insofar as large linewidths may be indicative of warm H I, and calls 
into question the likelihood of forming molecular clouds from colliding H I flows. We also find that the ratio of 
molecular to atomic gas shows no significant correlation (or anti-correlation) with the stellar surface density, 
though a correlation with midplane hydrostatic pressure P/, is found when the data are binned in P/,. The latter 
correlation largely reflects the increasing likelihood of CO detection at high H I column density. 

Subject headings: Magellanic Clouds — ISM; molecules — galaxies: ISM 



1. INTRODUCTION 

A key step in the cycling of gas between stars and the inter- 
stellar medium is the formation of molecular clouds, which in 
turn serve as the sites of further star formation. Understanding 
this process requires observations of both the atomic gas com- 
ponent, observed most readily using the 2 1 -cm line of H I, and 
the molecular gas component, observed most readily using 
the 2.6-mm 7=1— >0 line of CO. While the information pro- 
vided by these two spectral lines is not unambiguous — ^both, 
especially CO, can be affected by optical depth and excitation 
effects — they provide rough indications of gas column densi- 
ties, mean velocities, and velocity dispersions which can be 
used to constrain cloud formation models. 

Balancing H2 formation on grains and photodissociation 
by far-ultraviolet (FUV) radiation leads one to expect a core- 
envelope structure for molecular clouds, in which dense CO- 
emitting cores are shielded by H I envelopes. The ten- 
dency for H I mass surface densities to saturate at E ~ 
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10 Mq pc-2 or A^H ~ 10^' cm-2 (e.g., IWong & Blitj |2002|) 
suggests that much of the H I gas in galaxies may con- 
sist of such env elopes, with any add itional gas being trans- 
formed into H2 (IShava & Fede rman 1987). Direct evidence 
for H I envelopes ar oun d CO clouds has been p resented by 
IWannier etaH (Il983h and lAndersson et all (119911) . who found 
H I emission maxima near the edge s of CO clouds. More 
recently. Goldsmith & Li (20051 and [Goldsmith et al.l (l2007h 
have used time-dependent models of H2 formation, in con- 
junction with observations of H I self-absorption in molecular 
clouds, to show that compressed H I clouds become molec- 
ular from the inside out; the clouds they study have inferred 
minimum ages of 3-10 My r. Since the inte rstellar medium 
(ISM) is thermally bistable dField et al.l ll969). the raw mate- 
rial for molecular cloud formation is expected to be the cold 
neutral medium (CNM) traced by 21 -cm absorption studies, 
which occupies a much smaller volume than the low-density 
warm neutral medium (WNM). 

While molecular clouds are generally assumed to be quasi- 
equilibrium structures, a more dynamic view of the interstel- 
lar medium has recently been gaining prominence. In this 
picture, shocks driven by large-scale motions are the domi- 
nant process for producing density enhancements, and molec- 
ular clouds can form rapidly downstream of shock fronts (e.g., 
[Bergin et al. 2004; Glover & Mac Low 2007), a scenario con- 
sistent with the relatively short inferre d lifetimes for molecu - 
lar clouds in the solar neighborhood (iHartmann et al.ll2001h . 
If the timescale on which gas experiences shocks is shorter 
than the timescale for thermal instability to operate, then the 
CNM may not be a necessary precursor to molecular clouds, 
and one can even envision mol ecular clouds form ing directly 
from shocked WNM gas (e.g.. 'Bergin et al.'l 2004l). We note , 
however, that the simulations of Glover & M ac Low! (l2007l) . 
which form H2 on timescales of only ^1 Myr, begin with gas 
already at densities typical of the CNM (^100 cm"^). 

It is unclear how rapid molecular cloud formation would af- 
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feet the observed relationship between H I and CO. H I would 
still be a necessary ingredient for cloud formation, and in- 
deed some turbulent cloud formation models predict a long 
"latency" period of H I accumulat ion prior to a molecular 
cloud becoming observable in CO (IVazquez-Semadeni et alJ 
1200 7). Molecular clouds formed recently from colliding 
gas flows should be associated with larger H I line widths, 
but even this observational signature may be obscured if 
gravitational collapse is delaye d while H I flows are active 
(IVazquez-Sema deni et al.l 12007 ). or if the H I line width is 
dominated by the WNM, whose thermal velocity dispersion 
exceeds the tur bulent velocit y dispersion of the newly-formed 
CNM (IKovama & Inutsukal 12002) . Another possible signa- 
ture of turbulence may be the existence of a porous density 
structure with significant mixing of atomic and molecular ma- 
terial (e.g., Hennebelle et al. 2008), thus tending to erase any 
sharp column density threshold for H2 formation associated 
with self-shielding. 

The Large Magellanic Cloud (LMC) is the nearest star- 
forming galaxy to our own, and its nearly face-on aspect 
makes it an excellent target for studying the relationships be- 
tween atomic and molecular clouds and star formation. In 
this study we compare the global CO and H I distributions at 
2. '6 resolution, corresponding to a spatial scale of ~40 pc (as- 
suming a distance of 50 kpc to the LMC). While this scale 
resolves only the largest molecular clouds, it is comparable to 
the linear dimension of the area from which atomic gas must 
be gathered in order to form a giant molecular cloud of mass 
lO'^-lO^ Mq. We consider the detectability of CO emission 
and its dependence on H I column density, peak brightness, 
and velocity dispersion. We also evaluate the coincidence be- 
tween CO and H I emission in velocity space using a Gaussian 
decomposition method. Finally, we consider whether simple 
estimates of the hydrostatic disk pressure are able to predict 
the location of regions with high molecular to atomic gas ra- 
tios. 

Previous comparisons of CO and H I emission in the LMC 
have focused on specific regions or been undertaken at spa- 
tial scales much larger than 3'. The only global c omparison, 
asi de from a brie f discu ssion in lBlitz et al.l (120071) . was made 
by iMizuno et alj ( 12001b . who focused on the association of 
the large-scale "L" and "D" components of the galaxy seen in 
single-dish H I maps (15' beam FWHM) with the CO clouds 
mapped in the first NANTEN survey. In a companion paper 
(Fukui et al. 2009, hereafter Paper II), we will examine the ex- 
tent of the HI emission (in position and velocity space) around 
CO-detected clouds, and relate the CO and H I properties of 
giant molecular clouds with their star formation activity. 

2. OBSERVATIONAL DATA 

The H I data for the LMC were taken with the Australia 
Telescope Compact Array (ATCA) and the Parkes telescopes, 
both part of the Australia Telescope National Facility.'^ De- 
tails of the observations are given by Kim et al. ( 1998) and 
IStavelev^m ith et al. (2003), and the data combination is de- 
scribed by [Kim et al.. (2003). The resulting data cube covers 
a region of 7.°5 x 7.°5 at a spatial resolution of 1' and with a 
1(7 noise level of 2.4 K in a 1.65 km s"' channel. To facilitate 
comparison with the CO data, we regridded the H I spectra to 
the LSR velocity frame and smoothed the maps to a resolu- 
tion of 2'.6; the resulting RMS noise in the channel maps was 

The Australia Telescope is funded by the Commonwealth of Australia 
for operation as a National Facility managed by CSIRO. 



0.6 K per channel. At this resolution the maximum observed 
H I brightness temperature is ^110 K, so the peak signal-to- 
noise ratio is ^190. 

The CO data were taken with NANTEN, a 4-m radio tele- 
scope operated by Nagoya University at Las Campanas Ob- 
servatory in Chile. A region of ~30 square degrees was ob- 
served in 26,877 positions with a 2' grid spacing, slightly un- 
dersampling the telescope primary beam (HPBW=2f6). Each 
position was observed for about 3 minutes on-source, yield- 
ing a typical la RMS noise of 0.07 K in a 0.65 km s"' chan- 
nel. The maximum observed CO brightness was 2. 1 K, corre- 
sponding to a peak signal-to-noise ratio of ^30. A cube was 
formed from the spectra using the GLS projection in 1950.0 
coordinates from 150 to 350 km s"' in the vlsr frame. How- 
ever, spectra in some regions of the galaxy only cover 100 
km s"' in velocity, generally centered near the peak H I veloc- 
ity. Further deta ils about the obser vations and data processing 
can be found in Fukui et aTl ( l2008h . 

Since different parts of the LMC were observed to different 
sensitivity limits in CO, we constructed a noise image using 
the median absolute deviation (MAD), considering only chan- 
nels in each spectrum with T[, < 0.5 K. The MAD, which is 
more resistant to outliers than the standard deviation, is ob- 
tained by calculating the absolute values of the residuals (de- 
viations) from the median of the spectrum, then taking the 
median. The Icr noise estimate is then 1.48 x MAD. Obvious 
rectangular patches in the noise image indicate that large re- 
gions of the galaxy were observed to comparable sensitivity, 
presumably during continuous observing sessions. To empha- 
size these large scale noise variations (which span more than 
a factor of 2, from 0.045 to 0.1 K) and reduce the intrinsic 
scatter when deriving noise estimates from individual spectra, 
we applied a moving 3x3 pixel median filter to the noise 
image using the MFILTER task in GIPSY. 

The integrated H I and CO intensities were derived straight- 
forwardly by summing the intensity for all spectral channels 
at a given position. iFukui et alJ (l2008h found that all of the 
detected CO emission was restricted to a vlsr range of 210- 
310 km s"' . To reduce the noise contribution to the integrated 
CO map, we only summed the CO emission from 200-325 
kms"', although we summed the H I emission from 175- 
350 km s"' in order to accommodate the typically broader 
H I linewidths (ct,, w 12 kms"') compared to CO (ct,, ~ 3 
km s"'). Figure [T] shows the resulting integrated H I image. 
The thick solid contour represents the region over which CO 
emission has been observed. For most of our analysis we have 
avoided identifying and blanking emission-free regions, so as 
to provide a simple estimate of the sensitivity at each location. 
However, for purposes of overlaying a fiducial 1 K km s"' CO 
contour in Figure[T]we have used a signal detection algorithm 
(discussed in 33.3l below) to only sum channels which are ad- 
jacent to, or coincide with, 3(t detections. Visual comparison 
between the resulting CO contour and the distribution of pix- 
els detected in integrated intensity without blanking (shown 
as crosses) indicates that they are largely consistent. 

3. RESULTS 

3.1. CO Intensity vs. H I Properties 

In this section we consider the detectability of CO emission 
as a function of H I peak brightness temperature (hereafter 
Tmax), linewidth, and integrated intensity. We phrase the com- 
parison in this way because the limiting factor in comparing 
the two species is the ability to detect CO: not only is the peak 
signal-to-noise ratio a factor of ~6 higher for the H I data, but 
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Fig. 1 . — CO contours from NANTEN (at the 1 K km s ' level) overlaid on integrated H I emission smoothed to the resolution of the CO data. Units of the 
H I intensity are also K km s"' . To reject noise in the contour map, a blanking mask constructed with the CPROPS signal detection algorithm described in j|3.3l 
has been applied. The heavy solid contour represents the region observed with NANTEN. Small crosses indicate the pixels which are considered CO detections 
without the use of a blanking mask. The southeastern H I overdensity (SEHO) is identified at the lower left of the map. 



the fraction of spectra with T,^.^^ exceeding 50% of the global 
T^^„ is much higher for H I (8.8%) than for CO (0.15%), im- 
plying that extremely few regions in the galaxy are bright in 
CO. 

The comparison between CO and H I can be done in 2-D 
(using e.g., integrated intensity maps) or in 3-D (as a function 
of radial velocity as well as position). Each has advantages 
and disadvantages: a 2-D analysis can lead one to compare 
unrelated material as a result of line-of-sight superposition, 
whereas a 3-D analysis can weaken an underlying correlation 
if there are differences in intrinsic linewidths resulting from 
differences in temperature or optical depth. In this paper our 
main approach is a 2-D one, but in ^3.3l we employ multiple 
Gaussian fits to better isolate the relevant H I component. We 
defer a full 3-D analysis of the data cubes to Paper II. 

For the present analysis, we select every second pixel along 
each of the spatial axes, to ensure each pixel is statistically 
independent (pixel separation 4')- Figure |2] compares the 
integrated H I and CO intensities, with pixels detected in 



both H I and CO shown as filled circles, and pixels detected 
only in H I shown as upper limit symbols (the 3cr sensitiv- 
ity limit for detecting H I, /hi~33 K km s"', lies beyond the 
left edge of the panel). On the top and right plot axes we 
have converted the integrated intensities into equivalent hy- 
drogen nucleon column densities, assuming optically thin H I 
and a CO-to-Ha conversion factor appropriate to the Galaxy, 
MH,)//co=2 X 10^" cm-- (K km s"')-' iStrong & Mattoxl 
Il996l) . Upper limit symbols correspond to 3cr upper limits 
for the integrated CO intensity. Of approximately 6700 pixels 
in the region observed by NANTEN, only 215 (3.2%) are de- 
tected in both tracers, while only 22 (0.3%) are not detected in 
either tracer Only 1 pixel is detected in CO but not H I, cor- 
responding to the bright continuum source 30 Doradus where 
HI is seen in absorption. Thus, over most of the area surveyed 
by NANTEN, H I is detected but not CO. Below we discuss 
the detectability of CO as a function of several characteristic 
properties of the H I spectra. 

3.1.1. HI Integrated Intensity 
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Fig. 2. — Integrated CO intensity vs. H I intensity for 2' X 2' pixels separated by 4'. Black solid circles are CO detections, red arrows represent 3cr upper 
limits. At the bottom is a histogram of the CO detection fraction in equally spaced logarithmic bins of /hi, with a Hnear ordinate peaking at the value indicated. 
Since only about a quarter of the pixels are shown, dashed histograms indicate the range of values obtained from the other three sets of pixels. The diagonal solid 
line indicates a molecular mass fraction of 0.5, assuming optically thin H I and a Galactic value of the CO-to-H2 conversion factor. Correlation coefficients are 
explained in the text. 

but not sufficient for detecting CO. 

We find no evidence for a threshold H I column density 
above which the detection fraction suddenly increases, nor is 
there indication of a "saturation" H I column density above 
which additional gas appears purely in molecular form. In 
particular, there is a large range in logA^H (from 21-21.8) 
over which CO is detected, and the lack of CO detections 
at lower A^h may be a reflection of our sensitivity limit. As 
demonstrated below 03.4b . however, the azimuthally aver- 
aged value of Nh is close to 10^' cm"^ {^10 Mq pc"^) across 
a large range in radius, so one might easily conclude from 
the radial H I profile that this represents a "saturation" value, 
even though locally much higher H I column densities are ob- 
served. 

The error bars shown in the detection fraction histograms 
are derived as follows. First, the uncertainty in the number 
of pixels in a bin is assumed to arise from counting errors, 
so it is given by ^/N. The uncertainty in the number of CO 
detections in a bin is estimated as the number of pixels with 
values between 3.5cr and 4.5a, where a is the RMS noise of 
the CO map. Finally, the uncertainty in the detection frac- 
tion is given by standard error propagation assuming normally 



Figure |2] illustrates the tendency for strong CO emission 
(as measured by the integrated line intensity) to be associated 
with strong H I emission. The correlation is by no means per- 
fect: many regions with strong H I emission show little or no 
CO emission. The Spearman rank correlation coefficient, cal- 
culated for points detected in both tracers, is only rdet = 0.24, 
and the generalized rank correlation coefficient, taking into 
account CO upper limits using a survival analysis ( llsobe et all 
Il986l) . is only r^an = 0.12. Both correlations nonetheless ap- 
pear to be larger than would be obtained by chance, as each 
lies outside the range of values given in brackets, which repre- 
sents the full range of correlation coefficients obtained when 
randomly shuffling the abscissas of the data points 50 times. 

A histogram of the CO detection fraction is shown at the 
bottom of the plot. Placing the pixels into logarithmic bins 
of increasing H I intensity, the fraction of pixels with detected 
CO increases monotonically, peaking around 33% (although 
the last bin suffers from small number statistics). It is note- 
worthy that weak H I emission (below an integrated inten- 
sity /hi ~ 700 K km s"', or A^h ~ 10^' cm"^) is never asso- 
ciated with strong CO emission (above an integrated intensity 
^co ~ 3 K km s"'). Thus, high H I intensities are necessary 
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Fig. 3. — (a) Integrated CO intensity vs. H I peak briglitness temperature for 2' X 2' pixels separated by 4'. (b) Integrated CO intensity vs. H I velocity 
dispersion. For both panels, symbols and histograms are as in Figure|2] Values to the left of the vertical dot-dashed lines are below our adopted sensitivity limits 
and thus excluded. 
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FlG. 4. — (a) Integrated CO intensity vs. H I intensity for 8' X 8' pixels, (b) Integrated CO intensity vs. H I intensity for 16' X 16' pixels. For both panels, 
symbols and histograms are as in Figure|2] except that green solid circles indicate pixels which are classified as detections even though none of the constituent 2' 
pixels were considered detections. 



distributed errors. The dashed histograms represent the range 
of detection fractions observed when separately analyzing the 
four different pixel samples we have available, since we only 
select a quarter of the original pixels for analysis. The quoted 
errors appear to be consistent with these values, if not overly 
conservative. 

The slanted line in Figure |2] labeled /moi=0.5, indicates 
equality between atomic and molecular column densities (ex- 
pressed in terms of A^h)- Evidently the column density of 
molecular gas rarely exceeds the column density of atomic 
gas; the LMC appears to be an H I-dominated galaxy. (Al- 
though the upper limits on the left side of the figure formally 
allow high molecular fractions in this region of parameter 
space, we note that if CO emission actually existed at the 3(t 
level for all non-detected pixels, the total CO flux would ex- 



ceed the measured flux by a factor of 9). We caution that if 
the CO-to-H2 conversion factor is larger in the LMC than in 
the Galaxy, as expected from its low metallicity and dust con- 
tent, the molecular column densities will be underestimated, 
and the equality line will shift downward (e.g., by 0.6 dex for 
a factor of 4 change in A^h)- This would increase the number 
of molecule-dominated regions, though they would remain in 
the minority. On the other hand, the H I column density could 
also be underestimated if the H I opacity is sig nificant. We 
return to a discussion of these uncertainties in § ^4.1.1|j4.1.2l 

3.1.2. H I Peak Brightness Temperature 

High peak H I brightness temperatures also appear to be a 
requirement for CO detection, as indicated by Figure 13 a). A 
minimum cutoff level of rn]ax=2.4 K has been applied, cor- 
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responding to 3(7 in the H I cube. Bright CO is detected al- 
most exclusively in regions where T^nx > 20 K, and the like- 
lihood of detection increases monotonically with T^^ax- Once 
again, however, there is no indication of a threshold in T^^^ for 
CO detection, and a clear one-to-one correlation is lacking: 
many high-brightness H I filaments are not detected in CO. 
The Spearman rank correlation coefficient between rinax(H I) 
and Icq is only rdet = 0.15; including upper limits yields a sim- 
ilar value (rsm-v = 0.14), though with somewhat greater signif- 
icance. 

Although 7;„ax will scale with /hi for constant H I line width, 
it can also highlight regions of significant H I opacity: in the 
Galaxy, sight lines with Tmax > 60 K typically s how H I in ab- 
sorption with r ~ 1 jBraun & W alterbos 199^. iBraunI (Il997h 
has argued that Tinax images of nearby spiral galaxies serve 
as direct probes of the cool neutral medium (CNM), the H I 
phase most similar in temperature and density to molecular 
clouds. To check if beam dilution was responsible for the 
low values of Tmax observed towards some CO-detected re- 
gions, we extracted the 38 CO-detected pixels with T^ax < 30 
K (just 5% of all CO detections) and measured their H I peak 
brightness temperatures at the full 1' resolution of the ATCA 
data. The average value of T^ax at 1' resolution is 31 K, com- 
pared to 23 K at 2. '6 resolution. Thus, a small but significant 
number of regions appear to show detectable CO without high 
brightness H I, although the general conclusion remains that 
high values of T^^ax are necessary for detecting CO. 

3.1.3. HI Velocity Dispersion 

In contrast to the integrated intensity and peak brightness, 
the H I line width, as estimated by the 2nd moment of the 
line profile (excluding pixels below a 4(t threshold in /hi, and 
rejecting linewidths <3 kms"'), shows no relationship with 
CO [FigureOb)]. Thus, H I velocity dispersion is a poor pre- 
dictor of whether CO emission will be detected. One inter- 
pretation of this result is that colliding H I flows are unlikely 
to be the direct precursors to GMCs. Alternatively, the 2nd 
moment of the H I line profile may not reflect the local H I 
velocity dispersion at any position in space, due to blending 
of multiple gas components along the line of sight. This pos- 
sibility is addressed in 33.31 where we deduce a weak anti- 
correlation between cr,.(HI) and CO detectability following a 
Gaussian decomposition of the line profile. Our results indi- 
cate that if GMCs form from colliding H I flows, then much 
of the turbulent kinetic energy dissipates before CO becomes 
easily visible. 

3.2. Effect of Varying Spatial Resolution 

To determine whether the poor point-by-point correlations 
between CO and H I found in 33.1l result from a strong global 
correlation breaking down on small spatial scales, we spa- 
tially binned the original 2' pixels into 8' square and then 
16' square pixels and repeated the correlation analysis. Some 
improvement in the correlation is to be expected, simply by 
virtue of blending distinct emission peaks so as to empha- 
size the large-scale galactic structure. We indeed find larger 
correlation coefficients (rsuiv=0.25 and 0.37 respectively), but 
caution that the range of values spanned by the permutation 
tests also increases. We also see a modest increase in the frac- 
tion of the "total" (integrated over the entire datacube) CO 
flux recovered in the CO-detected pixels, from 54% to 67% 
(8' pixels) and 74% (16' pixels), although such measurements 
are sensitive to any errors in the spectral baseline. Finally, we 



note an increase in the fraction of CO-detected pixels (Fig- 
ureHli. This appears to reflect the basic tendency for detection 
within any of the pixels contributing to a bin to result in the 
larger, averaged pixel being considered a detection. In fact, 
of the 8' pixels observed in both CO and H I, only 0.4% were 
classified as detections when none of the 16 constituent 2' pix- 
els were so classified (these are shown as solid green circles in 
FigurelUl. Thus, the increased detection fraction after binning 
is not due to combining pixels that lie just below our detection 
limits. 

As indicated by Figure |4] measuring the correlation on a 
coarser grid still leads to a large range in Iqo at a given /hi, 
even as the overall CO detection fraction increases. In partic- 
ular, there is little evidence for a change in the tightness of the 
correlation above or below a particular averaging scale. This 
confirms the visual impression from Fig.[T]that the differences 
between the H I and CO distributions exist on large as well as 
small scales. Since most of additional flux that is recovered 
by spatial averaging appears to be located close to CO peaks 
already detected at high resolution, any additional flux tends 
to correlate with the H I emission in a similar way. 

3.3. Association in Velocity Space 

Our approach of integrating in velocity takes no account of 
whether H I and CO are actually coincident in velocity space. 
In Paper II, we will undertake a full 3-D comparison of the 
data cubes. However, we can test whether velocity integra- 
tion has strongly skewed our results by performing an au- 
tomated Gaussian decomposition of the H I cub e using the 
GAU FIT task in the MIRIAD software package (ISault et all 
Il995h . While such an analysis cannot guarantee that we are 
able to compare atomic and molecular gas at the same line- 
of-sight position, it makes it less likely that a correlation (or 
lack thereof) is solely due to projection effects. We fit up to 
two Gaussian components at each sky position, and if CO is 
detected at the same position, the H I component closer to the 
CO velocity (if within 20 km s"') is taken to be associated 
with CO. The component with the highest amplitude (peak 
Th) is taken to be the dominant H I component. For about 1/6 
of the pixels we find a significant improvement in when 
going from single to double Gaussian fits, especially in the 
central and southeastern parts of the LMC, w here two widely 
separ ated H I components have been noted (iLuks & Rohlfsl 
IT992h . 

For purposes of comparison with the H I Gaussian fits, we 
assume the CO profiles are single-peaked, which is a good as- 
sumption for the vast majority of spectra, for which fitting an 
additional component has little impact on x^. In order to min- 
imize noise contamination in determining the radial velocity 
at the peak of each CO spectrum, we applied a signal detec- 
tion algorithm adapted from the IDL-based CPROPS package 
of iRosolowsky & LeroyI (12006) . A masked datacube is gen- 
erated by starting at the 3a contour (required to encompass at 
least two consecutive channels and three adjacent pixels) and 
expanding it to include adjacent pixels in position and velocity 
above a la level. The CO velocity is then taken as the veloc- 
ity at the peak of each spectrum within the masked datacube. 
Of the 479 CO spectra in the full-resolution map that at least 
partially survive the masking process, we find that the vast 
majority (95%) associate closest in velocity with the domi- 
nant (or sole) H I component, with only 4% associated with 
the lower amplitude component. Only 1% of the CO detec- 
tions were unassociated with H I within 20 km s~' ; these spec- 
tra were usually of low signal-to-noise or found near 30 Do- 
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Fig. 5. — Scatter plots of the CO-H I velocity difference as a function of (a) CO intensity and (b) width of the associated Gaussian H I component. Open 
circles represent pixels where the dominant (highest peak 7^,) H I component is associated in velocity with the CO, while filled triangles represent pixels where 
the secondary component is associated. Contours indicate the density of points (including both components), and the lengths of vertical bars above the plotted 
points indicate the FWHM of the distribution. 
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Fig. 6. — Histograms of (a) H I intensity and (b) H I velocity dispersion for the entire region observed in CO. The top sub-panels show the distribution of all 
pixels using diagonal shading and the distribution of CO-detected pixels using black shading. The middle sub-panels show the same two distributions (as green 
and red shaded histograms respectively) but with the H I restricted to the dominant component in a two-Gaussian decomposition. The bottom sub-panels show, 
on a linear scale, the fraction of pixels detected in CO, with H I restricted to the dominant Gaussian component shown as solid red circles. 



radus, where the strong 1 .4 GHz continuum leads to strongly 
inverted H I profiles that cannot be fit by the Gaussian decom- 
position. 

Although by definition the velocity difference between as- 
sociated CO and H I components is less than 20 km s"', in 
most cases the difference is considerably smaller (Figure |5]l. 
In fact, for H I profiles best fit by two Gaussians, a sig- 
nificant fraction (19%) of the CO detections lie within 20 
km s"' of both H I Gaussian components. No clear relation 
between the velocity difference and the CO intensity is ap- 
parent [Fig.|3a)], with the distribution of velocity differences 
exhibiting a nearly constant width of ~10 kms"', indepen- 
dent of the CO intensity. Not surprisingly, many of the largest 
velocity differences occur when the associated H I component 



has a large line width, in which case the velocity centroid may 
be less clearly defined or several unresolved H I components 
may be superposed along the line of sight [Fig-IStb)]. 

How are the detection histograms in Figures |2] and [3] af- 
fected by the Gaussian decomposition? In the upper panels 
of Figure |6] we show the distributions of H I intensities and 
velocity dispersions over the full region observed by NAN- 
TEN and over just the CO-detected pixels. The middle pan- 
els show the equivalent distributions after isolating only the 
dominant Gaussian component (here the velocity dispersion 
is measured as the Gaussian width, which is equivalent to 
the 2nd moment). As expected, the Gaussian decomposition 
shifts both the H I intensity and velocity dispersion to lower 
values. In the bottom panels, a marginal increase in the CO 
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Fig. 7. — H I (top) and CO (bottom) channel maps integrated over three velocity ranges given with respect to a circular rotation model. Contour levels are 
400n K km s~' for H I and 2.4n K km s~' for CO, where ra=l,3,5,7. A cross indicates the adopted rotation center. 
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detection rate at the highest H I intensities is apparent follow- 
ing decomposition, underscoring the fact that CO associates 
with the dominant H I component. On the other hand, the CO 
detection fraction drops to zero for high H I velocity disper- 
sions, even in bins which remain highly populated. Thus, the 
Gaussian decomposition appears to reveal an anti-correlation 
between cr,. and CO detectability, although in all cases the CO 
detection fraction is low (<6%). The shift in CO detections 
toward lower cr,, suggests that many of the detections are in re- 
gions where the total line profile is artificially broadened (as 
measured by the 2nd moment) due to the presence of multiple 
H I components. 

To emphasize that CO emission primarily occurs near the 
velocity of the dominant H I component, regardless of how 
regular that component appears in velocity space, we have 
fit a circular rotation model to the H I line-of-sight veloc- 
ity field, estimated using the 1st moment of the line profile. 
Given the large angular size of the galaxy, we first correct for 
th e transverse motion of the L MC, following the prescription 
of Ivan der M arel et al. (2002), using the proper motion vec- 
tor derived by Kallivavalil et al. (2006) using Hubble Space 
Telescope imaging. We assume a disk inclination of 35° and, 
using the ROTCUR task in the GIPSY software package, de- 
rive a mean kinematic position angle of 340° and a kinematic 
center of 5'' 19"'5, -68°53^ (J 2000), about 14' from that de- 
termined by lKim et al.l ( |1998|) (this difference, while large, is 
within the formal errors of the fit). Shifting the H I spectra 
to the mean velocity of the circular rotation model yields a 
prominent low-velocity excess in the central and southeastern 
regions (Figure |7li, which can be identified as the "L" compo- 
nent of Luks & Rohlfs ( 1992). In addition, a (weaker) high- 
velocity component is also apparent, suggesting a bifurcation 
in the gas layer a t the e astern edge of the galaxy, as discussed 
by lNidever et al.l(l2008l) . Shifting the CO spectra according to 
the same model, we find that the CO is not exclusively asso- 
ciated with the main disk component: at the eastern edge of 
the galaxy, it is associated primarily w i th the L component, as 
was previously noted by Mizuno et al.l (1200 ll) . The CO excess 
at low velocities is also apparent in spatially averaged spectra, 
both before and after subtraction of the circular rotation model 
(Figure ISll. 

3.4. Molecular to Atomic Ratio 

IWong & Blit j (l2002h . iBlitz & Rosolowskvl (l200l l2006h 
and lLerov et al.l (l2008l) have examined H I and CO in nearby 
spiral galaxies and found that the ratio of molecular to atomic 
gas surface density correlates strongly with the hydrostatic 
pressure at the disk midplane. Pi,, estimated for a two- 
component disk of gas and stars. As given by Elmegreen 
(11989), Ph can be estimated from the observed mass surface 
densities of gas and stars (Egas and E*) as follows: 



ttG 
P;, = — E. 



5jgas"l E, 



(1) 



where the term in parentheses is an estimate of the total dy- 
namical mass within the gas layer, and ag and ct* are the ve- 
locity dispersions of gas and stars, respectively. 

To investigate the possible role of pressure in determin- 
ing the H2/H I ratio R,^oi in the LMC, we first consider two 
quantities contributing to Pj, that can be determined indepen- 
dently of the gas surface density. In Figure |9ja) the inferred 
molecular to atomic surface density ratio, calculated at pix- 
els where /hi > 110 K kms"', is shown as a function of 
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Fig. 8. — (a) Mean CO (histogram) and H I (curve) brightness temperature 
as a function of heliocentric velocity, calculated over the area of sky shown 
in Fig.|2] The CO data have been scaled up by a factor of lO'. (b) Same 
as (a), but after shifting the velocity scale of each spectrum to be centered 
on a circular rotation model. The dashed lines represent the boundaries of 
integration used for Fig. [7] 

the stellar surface density E», which dominates the gravita- 
tional potential in the inner disk. CO detections are shown 
as solid points, whereas the much larger number of CO upper 
limits have been binned and contoured. Thus, the contours 
should be thought of as the upper envelope of the CO non- 
detectio ns. We have used for E, the stellar surface density 
image of lYang et al.l (l2007i) . which is based on number counts 
of red giant and asymptotic giant branch stars as identified 
by 2MASS colors, and normalized to a total stellar mass of 
2x10'' Mq. No correlation between R^o\ and E» is apparent, 
nor is there any clear separation between the CO detections 
and non-detections. Similarly, as shown in Figure |9|b), only 
a weak anti-correlation is found when comparing R^o\ with 
the H I velocity dispersion, (Ti,(H I), calculated from the 2nd 
moment of the line profile. Using the Gaussian dispersion of 
the dominant H I component for ct,, does not significantly alter 
this result. 

In Figure |9jc) we correlate /?moi against /hi, on a pixel -by- 
pixel basis, reproducing Figure |2] but now with a substantial 
"tilt" that results from dividing by the abscissa. The detec- 
tions and non-detections are clearly separated, but no clear 
relationship is evident when taking the upper limits into ac- 
count. Evidently, the apparent anti-correlation seen in the de- 
tected points can be attributed to the difficulty in measuring 
/?moi when /hi is weak. Finally in Figure |9td) we correlate 
Rmoi against Pi,, with only H I considered in calculating Egas. 
This should be a reasonable approximation for the LMC be- 
cause the contribution of H2 to the total gas mass is relatively 
small unless an extremely large CO-to-H2 conversion factor 
is assumed. We plot radial mass surface density profiles for 
stars, atomic gas, and molecular gas in Figure [TO] to verify 
this; the inferred H2 surface density is generally a factor of 
5-10 less than the H I surface density. We assume a constant 
velocity dispersion of 9 km s"' for the H I, based on the aver- 



10 





Fig. 9. — Ratio of molecular to atomic gas surface density as a function of (a) stellar surface density; (b) H I velocity dispersion; (c) H I integrated intensity; 
(d) average midplane hydrostatic pressure, calculated as described in the text. The CO detections are shown as solid circles. Contours, spaced between 15% 
and 95% of the peak in increments of 20%, indicate the density of points for upper limits only, corresponding to CO non-detections. Correlation coefficients are 
calculated as for Figures|2j2] 



age dispersion of the dominant Gaussian component, and 20 
km s"' for the stars (van der Marel et al. 2002). Overall the 
results are similar to those shown in Figure |3c), with again 
an apparent anti-correlation likely driven by a strong tilt in 
the distribution of upper limits. We therefore find no clear ev- 
idence for a relationship between R^oX and either Sgas or P/,. 

An obvious limitation of this analysis is the preponderance 
of high upper limits on Rmo\- This reflects the comparatively 
poor sensitivity of the CO data, which precludes useful con- 
straints on the CO/H I ratio where H I emission is weak. By 
contrast, the adopted H I brightness cutoff when calculating 
Rmoh corresponding to a column density of A^h=2 x 10^** cm"^, 
only excludes 2% of the pixels in the region observed in CO. 
One could increase the signal-to-noise ratio on 7?,noi by az- 
im uthally averaging Sga s, as has been done for spiral galaxies 
by IWong & Blitj (l2002h . but using the radial profiles shown 



in Figure [TO] fails to reveal any meaningful trends between 
/?,noi and P/„ given that the large fluctuations in the radial CO 
profile are not reflected in either the H I or the stellar profile. 
In ^4.2| we average T^moi in bins of increasing P/, to provide an 
alternative comparison to the spiral galaxy results. 



4. DISCUSSION 



4.1. 



Interpreting the CO-H I Correlation 

We have found that while H I appears to be a necessary con- 
dition for the appearance of CO emission, the intensities of 
the two tracers are poorly correlated across the LMC. Given 
the expected core-envelope structure discussed in fl] this is 
somewhat surprising. In the absence of significant overlap be- 
tween clouds, an increased number of molecular cloud cores 
should be accompanied by an increased mass in H I envelopes. 
Conversely, high H I column densities should promote H2 for- 



Molecular and Atomic Gas in the LMC. 1. 



11 



radius (kpc) 
1 2 



c 
-a 

o 
to 

3 



O 




50 100 150 

radius (arcmin) 

Fig. 10. — Azimuthally averaged radial profiles of stellar. H I and H2 mass 
surface density for the LMC. A Galactic value for the CO-to-H2 conversion 
factor has been assumed. The adopted center is 5''19l"5, -68°53' (J2000) 
and the adopted inclination and major axis position angle are 35° and 340° 
respectively. 

mation because more gas can in principle be shielded from 
external UV radiation. The absence of a strong correlation 
leads to three general classes of interpretation, each of which 
can be linked with one or more physical mechanisms. 

The first possibility is that the scatter in the correlation 
results from a partial depletion of H I as it is converted to 
H2. This could generate a large observed range in the ratio 
of CO to H I intensi ty, depending o n how far the conversion 
proceeds. Recently, iOtt et alJ (12008 ) have reported evidence 
of H I to H2 conversion, in the form of offsets between CO 
and H I peaks, for molecular complexes south of 30 Doradus 
mapped at 45" resolution. However, the larger scales probed 
by our analysis (by a factor of 4 in each linear dimension) 
should dilute the effects of such localized conversions, be- 
cause the H I appears to dominate the gas mass both locally 
(Figure |2]l and on large scales (Figure [Toll. Formation of H2 
therefore seems unlikely to have had a major effect on the H I 
intensity at the scales we are probing. 

The second possibility is that our assumption that CO and 
H I trace molecular and atomic column density faithfully is 
flawed. For instance, the dense gas may be incompletely sam- 
pled by CO emission, due to the time required for CO to be- 
come visible or variations in the degree of shielding against 
dissociating FUV radiation, or may be underluminous in the 
H I line due to high optical depth. We discuss these possibili- 
ties in greater detail in § ^4.1. 1 I and l4. 1 .2l below. 

The third possibility is that volume density and column den- 
sity are decoupled: large H I column densities do not necessar- 
ily imply the large volume densities characteristic of molec- 
ular clouds. This decoupling would arise naturally from a 
change in the scale height of the gas layer — due to fluctua- 
tions in the amplitude of turbulent motions or a variation in 
the gravitational potential — or may reflect differences in den- 
sity contrast which are smoothed over when averaging over 
^--^40 pc scales. We discuss these possibilities in greater detail 



in § ^4T3] and l4T4l below. 

4.1.1. Variations in the CO-to-Hi Ratio 

The use of CO as a tracer of H2 is usually justified with 
the following argument: while the CO(l-O) line is optically 
thick, the integrated line emission arises from an ensemble of 
virialized clouds that can be characterized by a common value 
of excitation temperature and density, and t hat do not shadow 
each other in position and veloc ity space (iMalonev & Blaclj 
[1988; Papadopo ulos et aT] l2002h . Even adopting these as- 
sumptions, however, it is unclear how reliably CO intensi- 
ties can be translated into H2 column densities for individual 
clouds, which would be expected to show variations in tem- 
perature, density, and degree of virialization. The role of op- 
tical depth variations is also unclear For clouds which are 
less optically thick, the observed emission may come from 
deeper, warmer layers, actually increasing the CO bright- 
ness. While we lack the data to address these issues prop- 
erly, we note that multi-transition studies of LMC clouds have 
generally found line intensity ratios t hat are consistent with 

CO having sign i ficant opacity (e.g., [Johansson et al .l 1 19941; 
[Israel et al.ll2003h . [Fukui et al.1 ([20081) find a strong correla- 
tion between CO luminosity and virial mass for the NANTEN 
cloud sample, and thus derive a CO-to-H2 conversion factor 
of 7 ± 2 X 10^° cm"^ (K km s~')"', roughly twice the value 
obtaine d from application o f the virial theorem to Galactic 
clouds ([Solomon et al.[[T987l) . and more th an three times the 
value obtained from Galactic 7-ray studies ([Strong & Mattox[ 
[T996h . 

A particular concern with regard to low-metallicity systems 
such as the LMC is that the CO-emitting region of a molecular 
cloud could become significantly smaller than the cloud itself, 
despite the high opacity of the CO line, as a result of CO pho- 
todissociation. Since dissociation of H2 occurs primarily at 
specific UV wavelengths (the Lyman- Werner bands), H2 can 
shield itself from photodissociation at sufficiently high col- 
umn densities (>10'^ cm"^). Much higher equivalent column 
densities are required for CO self-shielding to be important, 
however, due to its much lower abundance, making shielding 
by dust important for the survival of CO. Calculations there- 
fore suggest that CO may not be the primary re servoir of car- 
bon in some molecular regions (Ivan Dishoeck & Bl ack 1988). 
Especially in regions permeated by strong UV radiation fields, 
H2 may be present without abundant CO, weakening any 
correlation between H I and CO emission that might other- 
wise exist. Maloney & Black (1988) and Paketal. (1998|) 
have presented models showing how the CO-emitting parts 
of molecular clouds can be embedded in much larger H2 en- 
velopes. The importance of this effect is a strong function of 
the dust abundance (and thus metallicity) as well as radiation 
field. 

An examination of the location of bright H I not associated 
with CO emission (Figure [TTT i lends some support to this in- 
terpretation. Contours delineating regions with 7'max(H I)>60 
K generally intersect or enclose CO-detected pixels; this im- 
plies that bright H I emission is generally found adjacent to, 
if not coincident with, CO emission. Thus, many regions of 
strong H I emission may harbor H2 clouds but be insufficiently 
shielded from dissociating FUV radiation to permit survival of 
CO. The strong concentration of HI near 30 Doradus may be 
an example of this effect. Interestingly, the largest high-7;„ax 
complexes appear to be in the outer regions of the galaxy. An- 
other contributing factor may be differences in the formation 
timescales of H2 and CO, which can create variations in the 
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Fig. 1 1. — Image of peak H I brightness temperature, with regions where CO is detected (as determined by the CPROPS masJf described in j|3.3t excluded 
(shown as white regions). Contours are spaced by 20 K starting at 60 K. The location of the 30 Doradus star-forming region is indicated by a star, and the 
kinematic center is indicated by a cross. 



CO/H2 ratio across a cloud. Models of molecular cloud for- 
mation behind shock waves indicate timescales of ^10'' yr 
for CO to become abunda nt, several times lo nger than the for- 
mation timescale for H2 (Ber gin et al.ll2004) . Again, this is a 
result of the more stringent shielding requirements for CO for- 
mation. However, this effect ought to be suppressed when av- 
eraging over areas much larger than an individual cloud (^^10 
pc), as we have done in this study. 

Observations of H2 in absorption provide the most direct 
constraints on the amount of molecular material that is missed 
by CO surveys. Unfortunately, since such observations are 
performed at ultraviolet wavelengths, they are limited by the 
number of suitable background sources and a bias towards 
lines of sight with relatively low extinction. Analysis of the 
existing FUSE data ( Welty et al., in preparation; see also 
iTumlinson et alJl2002l) yields relatively low H2 column den- 
sities (A^H ^ 10^' cm"^) compared to those traceable by CO 



emission, and thus the low CO(l-O) detection rates found to- 
ward those sight lines are to be expected. Although there is 
always ambiguity in comparing emission and absorption mea- 
surements because they generally probe different volumes and 
different line-of-sight depths, deeper CO integrations towards 
these positions would be worthwhile. In addition, observa- 
tions of H2 (or suitable proxies) in absorption towards more 
highly reddened sight-lines near CO clouds would allow a 
more definitive comparison between absorption and emission 
studies. 

Recently, indirect evidence for an H2 compo nent outside 
of CO -emitting regions has been pres e nted i 7y iLerov et al] 
( l2007h for the SMC and [Bernard etal] ( l2008h for the LMC 
based on far-infrared measurements with Spitzer, which 
would otherwise imply a discrepancy in dust abundance be- 
tween low and high column density regions . On the other 
hand, UV absorption line measurements (ICartledge et alJ 
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I2005L and references therein) have inferred just such a dis- 
crepancy by comparing E(B-V)/N(\{ I) and £'(B-y)/A^(H2) 
toward multiple UV sight lines. Indeed, one might expect a 
higher dust abundance in H2 -dominated regions because of 
the important role of dust grains in H2 formation. Thus, it 
remains unclear to what extent an unseen H2 component is 
required by observational data. 

4. 1 .2. Variations in H I Optical Depth 

For an optically thick line profile, the integrated intensity 
will underestimate the column density. In principle, high 
opacity in the H I line could account for the different CO in- 
tensities detected towards regions of similar integrated H I in- 
tensity. Several studies of H I emission and absorption along 
sight lines pa ssing throu gh the LMC have been conducted 
jDickey et al.l 17994; Mebold et al.l ll997l: iMarx-Zimmer et all 
I2OOO), indicating widespread cool H I gas with 7=30-40 K 
and peak optical depths of r=0. 5-1.5. In par ticular, three 
of the lines of sight reported by iDickev et al.l (|T994) show 
bright CO emission: they ai-e towards 0521-699, 0539-697, 
and 0540-697, and show peak H I optical depths of 0.91, 
1.68, and 1.39 respectively. Since the H I column den- 
sity i s underestimated by a factor of roughly r/[l -exp(-r)] 
(e.g.. lDickev et al.ll2003h . towards optically thick sight lines it 
should be corrected upward by a factor of 1.3-2. This is in fact 
an overestimate because opacity in the line wi ngs is much less 
than at the line center (see e.g.. Fig. 2 of M arx-Zimmer et all 
12000) . Corrections of this magnitude would not significantly 
affect the trends seen in Figure |2] 

We also note that peak values of Th{l{ I) reach values of 
Tmax ^ 100 K (Fig.O; cool H I at temperatures of ~30 K (as 
inferred from previous absorption studies) would be a minor- 
ity constituent of the total H I brightness temperature, and thus 
not a major factor in the 2-D analysis considered here. 

4.1.3. Variations in Gas Thickness 

The poor correlation between CO and H I intensities could 
arise from a decoupling of volume and surface densities as a 
result of variations in the gaseous scale height. The thickness 
of the gas layer in a galaxy is related to both the pressure sup- 
port of the gas and the gravitational force acting on it. An 
increase in the gas scale height is most easily achieved by 
either increasing the velocity dispersion of the gas or decreas- 
ing the gravitational force on it. Consequently, if a thicker 
gas disk were responsible for lowering the volume density in 
H I-rich, CO-poor regions, one would expect these regions 
to exhibit higher H I velocity dispersion and/or lower stellar 
surface density. Fig. |9|a), however, shows no tendency for 
low H2/H I ratios to be associated with weak stellar gravity 
(low E*). Fig. |9|b) shows a hint of an anti-correlation be- 
tween /?moi and the H I velocity dispersion, although the trend 
is weak. If real, such an anti-correlation would be consistent 
with a lowering of the volume density due to increased tur- 
bulent support of the H I layer. Consistent with this idea, we 
noted in Figure|6|b) an apparent decrease in CO detectability 
at larger values of a^iVi I). 

On the other hand, turbulence would be expected to si- 
multaneously increase local volume densities by generating 
small-scale density structure (as discussed further in ^4.1.41 
below), possibly leading to increased CO emission. An en- 
ergy source for the turbulence that thickens the gas disk is 
also needed, which c ould be either ma ssive star formation 
or galactic dynamics dElmegreen & Scalo .2004.) . The abun- 
dance of supershells and the lack of strong density wave struc- 



tures in the LMC suggests that star formation should domi- 
nate, but star formation is expected to correlate with higher, 
not lower, gas density. Adding to the ambiguity are the pos- 
sible contributions of cosmic ray, magnetic field, and exter- 
nal pressures. These factors may not be negligible, especially 
given the interaction of the LMC with the Galactic halo, and 
their influence on the gas layer thickness is poorly understood. 

4.1.4. Variations in Density Contrast 

Two regions of similar average H I column density could 
have very different density structures when observed at higher 
resolution. A high degree of clumping, leading to larger den- 
sity contrasts, could increase the amount of material at high 
enough densities to radiate CO line emission, while still main- 
taining a fixed column density when averaged over a telescope 
beam. Conversely, a more uniform distribution of material 
could lead to similar average column densities but with rela- 
tively little high-density gas. 

A straightforward way to model density contrasts is to con- 
sider the probability distribution function (PDF) of the vol- 
ume density, which in the presence of turbulent forcing (e.g., 
IWada & Normanir2001h approaches a log-normal function. 
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IWada & NormanI (|2007|) show that the width a of the PDF 
depends on the ratio of the volume averaged density (p) to 
the characteristic density po, 
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They therefore derive the mass fraction of gas above a critical 
density pc to be 



fc=\i 1-erf 



ln(p7(p))-aV2 
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which increases strongly with <t. Adopting an average gas 
density (p) ^ 0.05 Mq pc~^ or nn^2 cm"-' for the LMC (us- 
ing a surface density of Sh ~ 10 Mq pc"^ and disk thickness 
of 200 pc, e.g. Kim et al. 1999), roughly 17% of the gas mass 
will have densities exceeding 100 cm"^ for a=2, but only 3% 
for (7=1.5. Therefore, if molecular clouds form primarily as a 
result of turbulent compression, their existence could be deter- 
mined primarily by variations in a, and high column density 
would be necessary but not sufficient for high volume density. 
For isothermal turbulence, the wi dth a of the PDF is a sim ple 
function of the rms Mach number dPadoan & Nordlundl2002) . 
althou gh such a relation may b e inapplicable to the ISM as a 
whole lWada & Normanll2007h . 

On the other hand, attributing changes in the fraction of 
dense gas to the effects of turbulence would lead one to ex- 
pect larger H I line widths in CO-detected regions, contrary 
to our results 03.1.3I ). A more plausible explanation is a sig- 
nificant component of warm, diffuse gas, which forms a sec- 
ondary peak di stinct from the lo g-normal portion of the den- 
sity PDF (Wada & Normanll200ll 12007). Despite its low den- 
sity, the WNM can make a large contribution to the column 
density averaged over a beam by virtue of its high volum e 
filling factor: in the simulations of IWada & NormanI (l2007h . 
80-90% of the volume is occupied by low-density gas. In 
addition, WNM profiles will be broader because of the much 
higher temperature of the WNM (factor of ^100 higher than 
CNM, corresponding to a factor of ^10 increase in thermal 
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Fig. 12. — Ratio of molecular to atomic gas surface density as a function 
of midplane hydrostatic pressure, averaged over all pixels in logarithmically 
spaced bins of P;, . The open triangles correspond to the southeastern region 
(SEHO) outUned in Figure[T] while the filled circles (to which a power law of 
slope 0.86 has been fit, as shown by the dotted line) correspond to the rest of 
the galaxy. Note that the highest pressures are only found in the SEHO. Error 
bars represent the standard error in the mean for each bin, and non-detections 
have been included in the bin averages. 

linewidth). It is therefore plausible that regions of relatively 
high H 1 column density but lacking in CO emission have a 
higher proportion of WNM gas, which would also account 
for the observed anti-correlation between CO intensity and 
H I velocity dispersion. Of course the physical mechanisms 
governing the CNM/WNM phase balance would still need to 
be identified; these would have to be more complex than a 
simple dependence on HI column density. 

Changes in density contrast in the ISM should lead to ob- 
servable changes in H I line brightness when observing at finer 
spatial resolution. If atomic and molecular gas are part of 
a continuum in gas density, one might also detect analogous 
signatures by comparing CO with high-density molecular line 
tracers such as HCO"^ and HCN, since the width of the density 
PDF will affect that fraction of material able to emit strongl y 
in those lines (e.g.. iBallesteros-Paredes & Mac Lowll2002l) . 
although abundances and interstellar chemistry will also af- 
fect their strength. Future observations should be able to test 
these possibilities. 

4.2. Molecular to Atomic Ratio and Pressure 

Equilibrium analysis of H2 self-shielding suggests that the 
column density of the H I envelope surrounding a molecu- 
lar cloud should scale approximately as 



Nn-cc 



urf 



-1.4 



(4) 



where «gas is the mean gas density, Ru is the rate constant for 
H2 formation on grains in units of cm^ s"' (depending pri- 
marily on metallicity), and Gsurt is the H2 photodissociation 
rate at the cloud surface ( Federman et al..,1979) . As the am- 
bient pressure (and thus cloud density) increases in the cen- 



tral regions of a galaxy, due to the deeper gravitational po- 
tential, molecular cloud enve lopes become thi nner and the 
H2/H I ratio (Rmoi) increases dElmegreenl 1 1 9891) . This trend 
has be en confirmed in spi ral galaxies with strong CO emis- 
sion ( Wong & Blitz 2002; Blitz & Rosolowsky 2006). Em- 
pirically, the correlation between hydrostatic pressure and the 
CO/H I ratio arises because the latter is a strong function of 
radius and also correlates well with the variation in the stellar 
surface density. However, in the case of the LMC, many of 
the massive CO complexes are found well outside the region 
of high stellar density (the bar). As indicated by Figure |9l a), 
no significant correlation between R^oi and is apparent. 

Figure |9id), on the other hand, does show an apparent sep- 
aration between detections and upper limits when comparing 
Rmoi and P/,. To investigate this further, we plot in Figure [12] 
the average values of R^o\ in logarithmically spaced bins of 
Pi,. Both detections and non-detections have been included 
in the averaging, and we have separated the southeastern re- 
gion of bright H I [dubbed the SE H I overdensity (SEHO) by 
iNidever et al.l (|2008)], as defined in Figure [T] from the rest of 
the galaxy. Excluding the SEHO, we find a good correlation 
between /J^oi and Pi,, with a best-fit power-law slope of 0.86. 
The slope is similar to that f ound bv lWong & Blitz! ilOOl ) and 
iBlitz & Rosolowskvl (l2006h . although the agreement may be 
misleading since both earlier studies only considered CO de- 
tections. Indeed, the correlation can largely be attributed to 
the increasing CO detection rate with increasing /hi seen in 
Figure |2] since most of the variation in Pi, results from vari- 
ation in Ehi [compare Fig.|9lc) with Fig.|9ld)]. Thus, while 
there is clearly a relationship between /?nioi and P/„ we pre- 
fer to attribute it to the fact that significant H I intensity is a 
prerequisite for CO emission, rather than any underlying re- 
lationship involving the stellar potential or turbulent gas pres- 
sure. 

The additional implication of Fig. [12] that /^moi may be- 
have differently in the SEHO region, is noteworthy given that 
this region is likely experiencing external ram pressure due 
to the LMC's motion in the Galactic halo, w hich may be 
compressing the gas disk jde Boer et al.lll99^ . X-ray spec- 
tra in this region suggest thermal pressures of ~ 10^ cm"^ K 
(Blondiau et al. 1997). For comparison, equation ([U gives a 
fiducial estimate for Pi, in the central regions of the LMC (us- 
ing the surface densities in Figure [TOli of P /k ^ lO'* cm"^ K, 
comparable to estimates for the solar neighborhood. Clearly, 
further examination of this and other interacting galaxy envi- 
ronments is needed to better understand how R^oX is governed. 

Differences in the relationship between Pi, and R,^o\ in the 
LMC as compared to molecule-rich spirals probably indicate 
that Pi, responds more sensitively to atomic gas in the former 
case and molecular gas in the latter case. The expression for 
Pi, given in equation ([TJ can be written as 



where 



Q* = 



Ph — T^gasCj 



2Q. 



(5) 
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and geff = 



gas 
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Here Qeff is an approximate gravitational instability parame- 
ter for a combined disk of gas and stars (Wang & Silk 1994). 
Note that Qeff/K is a measure of the disk dynamical timescale 
(e.g., Krumholz & McKee 2005), soP,, = T,„.^,ag/(2tay„). (The 
equivalent expression Pi, = pgcrl then yields a gas scale 
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height of hg ^ T,gns/Pg ^ crgkyn, as expected.) In the in- 
ner disks of CO-bright spiral galaxies, Sgas is dominated by 
molecular gas, a-g is relatively constant, and the dynamical 
timescale Q/k is expected on theoretical gr ounds to be pro- 
portional to the star formation timescale ( ElmegreenI l2002t 
iKrumholz & McK ee 2005), which for these regions is also 
fairly constant ([Wong & Blit3l2002l: llerov et alj |2008). Thus 
f/i scales linearly with the H2 surface density and with R^o\ 
(give n the saturation of the H I column density near 10 Mq 
pc"-, Lerov et aLll2008l) . In the case of the LMC, however, 
most of the gas mass is in H I, and furthermore the H I column 
density is generally not saturated, especially when measured 
locally. These conditions will ensure that f;, correlates pri- 
marily with Ehi, which results in a much weaker relationship 
with/?moi- 

Finally, attempts to predict R,^oi based on equilibrium val- 
ues of surface density and/or pressure may be frustrated if, 
as may be the case in molecule-poor environments, the time 
to reach equilibrium H2 abundances exceeds the lifetime of 
clouds against dissociation. Such a situation has been sug- 
gested by Leroy et al. (2007) in relation to the SMC. Under 
such conditions no equilibrium would be obtained, and the 
formation of H2 may be primarily dependent on local con- 
ditions such as dust abundance. Future comparisons of the 
H2 abundance in the LMC, as measured by UV absorption 
towards background sources, with other diagnostics of inter- 
stellar conditions could be used to test this hypothesis. 

5. SUMMARY AND CONCLUSIONS 

In this paper we have examined correlations between CO 
integrated intensity and H I integrated intensity, peak bright- 
ness temperature, and linewidth across the LMC. Our main 
results can be summarized below: 

1 . Large H I integrated intensities and peak brightnesses 
are prerequisites for detection of CO, with the detection 
fraction rising steadily for larger values of both quanti- 
ties, though never reaching 100% on spatial scales of 
-40 pc. 

2. On the other hand, neither condition is sufficient to 
guarantee CO detection, and the quantitative correla- 
tions of both quantities with CO intensity are poor. 

3. The large range in CO intensities at the highest H I in- 
tensities persists even when the data are spatially av- 
eraged on >200 pc scales, suggesting that local con- 
version of atomic to molecular gas is not the principal 
reason for the scatter 

4. The H I linewidth shows no clear correlation with CO 
intensity, suggesting that molecular clouds do not form 
directly from the collision of H I clouds. Decomposing 
the H I profiles into up to two Gaussian components and 
identifying the component more closely associated in 
velocity with CO actually suggests an anti-correlation 
between H I linewidth and CO intensity. This may pose 
a problem for turbulent models for forming molecular 
clouds, although it may also reflect a time delay be- 
tween the initial cloud compression and the appearance 
of widespread CO emission. 

5. While our analysis does not investigate possible differ- 
ences between CO and H I spectral profiles, we find 
from our simple Gaussian decomposition that the CO is 



almost invariably associated with the highest brightness 
H I component, irrespective of whether this component 
exhibits "disk-like" or "anomalous" kinematics. A full 
3-D analysis will be presented in Paper II. 



6. We find no clear correlation between the molecular to 
atomic gas ratio, as traced by /co//hi, and the hydro- 
static pressure Pi,, when examined on a point-by-point 
basis. Since this comparison is problematic because 
both quantities involve /hi, we also search for corre- 
lations between Ico/Im and the stellar surface density 
and gas velocity dispersion, but again find no signifi- 
cant trends. That the LMC does not follow the trends 
observed in CO-bright spiral galaxies may be a conse- 
quence of its relatively small molecular fraction. 

The fact that CO is always associated with significant H I 
emission, while not surprising, confirms the importance of 
atomic gas as the raw material for molecular cloud forma- 
tion. That strong H I emission is not always associated with 
CO suggests several possibilities, such as extensive H2 gas not 
traced by CO, variations in the timescale for CO formation, or 
a decoupling between surface and volume density as a result 
of changes in the gas layer thickness or density contrast. Al- 
though none of these possibilities can be ruled out entirely, 
we suggest that the presence of a low-density WNM which 
makes a significant contribution to the observed column den- 
sity along some lines of sight provides a natural explanation 
for regions exhibiting similar A^(H I) but different amounts 
of CO emission. In this context, the broadening of the H I 
linewidth in WNM-dominated regions might account for the 
weak anti-correlation between linewidth and CO emission. 

A widespread component of low-density H I would also 
have implications for predicting the molecular to atomic gas 
ratio, since such a component does no t contribute to shield- 
ing H2 against photodissociation (see IKrumholz et al.l l2008l 
for further discussion). Models which attempt to explain the 
balance between H2 and H I in terms of shielding would need 
to take this component into account. In addition, the time is 
ripe for undertaking a joint analysis of the existing H I emis- 
sion and absorption data in order to model the large-scale dis- 
tributions of the WNM and CNM in the LMC. The simple 
Gaussian fitting presented in this paper is inadequate for this 
task, given the complexity of the H I spectral line profiles, 
and a combination of spatial a nd spectral decomposition, as 
has been recently presented by iKim et aP (l2007h . appears to 
be a more promising approach. 

Finally, the utility of CO as a tracer of H2 should con- 
tinue to be investigated with improved data, especially for 
the M agellanic Clouds. The ongoing Mopra survey (Ott et aD 
l2008h will map CO across the LMC at an angular resolution 
of ~45", surpassing that of the H I synthesis data. Further 
constraints on the amount of H2 not associated with CO is 
now being provided by maps of far-infrared or sub-millimeter 
dust emiss ion, or dust extinction in the near- to mid-infrared 
(e.g., I mara & Blitzll2007HDobashi et al.ll2008i) . We have also 
emphasized the complementary role played by UV absorp- 
tion studies. In addition, further constraints on the physical 
properties of the CO-emitting clouds, through observations 
of high-excitation transitions or optically thin isotopologues, 
will help to clarify the conditions under which CO is easily 
detected. 
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